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ABSTRACT 

We perform a definitive analysis of Bianchi Yllh cosmologies with WMAP observations of 
the cosmic microwave background (CMB) temperature anisotropies. Bayesian analysis tech- 
niques are developed to study anisotropic cosmologies using full-sky and partial-sky, masked 
CMB temperature data. We apply these techniques to analyse the full-sky internal linear com- 
bination (ILC) map and a partial-sky, masked W-band map of WMAP 9-year observations. 
In addition to the physically motivated Bianchi Yllh model, we examine phenomenological 
models considered in previous studies, in which the Bianchi YIl h parameters are decoupled 
from the standard cosmological parameters. In the two phenomenological models considered, 
Bayes factors of 1.7 and 1.1 units of log-evidence favouring a Bianchi component are found 
in full-sky ILC data. The corresponding best-fit Bianchi maps recovered are similar for both 
phenomenological models and are very close to those found in previous studies using earlier 
WMAP data releases. However, no evidence for a phenomenological Bianchi component is 
found in the partial-sky W-band data. In the physical Bianchi Yllh model we find no evidence 
for a Bianchi component: WMAP data thus do not favour Bianchi Yllh cosmologies over 
the standard A Cold Dark Matter (ACDM) cosmology. It is not possible to discount Bianchi 
Yllh cosmologies in favour of ACDM completely, but we are able to constrain the vorticity of 
physical Bianchi Yllh cosmologies at (oj/H) < 8.6 x 10" 10 with 95% confidence. 

Key words: cosmology: cosmic background radiation - cosmology: observations - methods: 
data analysis - methods: statistical. 



1 INTRODUCTION 

The A Cold Dark Matter (ACDM) cosmological concordance 
model has recently emerged as an accurate description of our Uni- 
verse. In this model the current Universe is dominated by dark 
energy and dark matter, with structure seeded by primordial den- 
sity perturbations generated during an inflationary phase in the 
early Universe. Support for ACDM is derived from a range of re- 
cent high-precision cosmological observations, with measurements 
of the cosmic microwave background (CMB), in particular those 
made by the Wilkinson Microwave Anisotropy Probe (WMAP) 
(Komatsu et al. 2011; Larson et al. 2011; Hinshaw et al. 2012), 
playing a leading role. While the concordance model is undeniably 
successful, it is important to test the assumptions on which it is 
based. One of the most fundamental assumptions of ACDM cos- 
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mology is the cosmological principle; namely, that the Universe is 
homogenous and isotropic on large scales. Evidence that this as- 
sumption is inaccurate would necessitate revision of ACDM. 

In this article we seek to test the global isotropy of the Uni- 
verse. When studying phenomena beyond cosmological concor- 
dance it is important to do so in the context of a well-motivated 
cosmological model. We focus on the homogenous but anisotropic 
Bianchi models. In these models the assumption of isotropy about 
each point in the Universe is relaxed, yielding more general so- 
lutions to Einstein's field equations. For small anisotropy, as de- 
manded by current observations, linear perturbation about the stan- 
dard Friedmann-Robertson- Walker (FRW) metric may be applied, 
leading to a subdominant, deterministic contribution to the CMB 
fluctuations. In this setting CMB fluctuations may be viewed as the 
sum of a deterministic Bianchi contribution and the usual stochastic 
contribution that arises in the ACDM model. 

The induced CMB temperature fluctuations that result in the 
homogenous Bianchi models were first studied by Collins & Hawk- 
ing (1973) and Barrow et al. (1985), ignoring the effects of dark 
energy. Barrow et al. (1985) focus on the most general Bianchi 
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VII/j and IX types, corresponding to open/flat and closed universes 
respectively. Since tight constraints have already been placed on 
Bianchi IX models by Barrow et al. (1985), we follow the recent 
trend in the literature and focus on the open/flat Bianchi VII/j mod- 
els, in which geodesic focusing produces spiral-type contributions 
in the CMB. Although we focus on CMB temperature fluctuations 
only in this article, we note that the CMB polarisation contribu- 
tions induced in Bianchi models have also been studied recently 
(Pontzen & Challinor 2007; Pontzen 2009; Pontzen & Challinor 
2011). 

Bianchi Vll h models were first compared to COBE data by 
Bunn et al. (1996) and Kogut et al. (1997), and to WMAP data 
by Jaffe et al. (2005, 2006b). A statistically significant correlation 
between one of the Bianchi Vllh models and the WMAP internal 
linear combination (ILC) map (Bennett et al. 2003) was discov- 
ered by Jaffe et al. (2005) by modelling the CMB as an unknown 
Bianchi component on top of the (fixed) best-fit WMAP ACDM 
cosmology. However, it was noted that the parameters of the best- 
fit Bianchi component were incompatible with those of ACDM. 
Nevertheless, quite remarkably it was found that when the WMAP 
data were 'corrected' for the best-fit Bianchi map, some of the so- 
called 'anomalies' reported in WMAP data disappeared (Jaffe et al. 
2005, 2006b; Cayon et al. 2006; McEwen et al. 2006). A modi- 
fied template-fitting technique was performed by Land & Magueijo 
(2006) and, although a statistically significant template fit was not 
reported, the corresponding 'corrected' WMAP data were again 
free of many large scale 'anomalies'. 

Following this renewed interest in Bianchi models, the CMB 
temperature fluctuations induced in Bianchi Vll h models incorpo- 
rating dark energy were derived (Jaffe et al. 2006c; Bridges et al. 
2007). In this scenario, a degeneracy between the matter and dark 
energy densities, Q m and £l A respectively, is introduced, but the 
cosmological parameters of the best-fit Bianchi template found 
in WMAP data nevertheless remain inconsistent with constraints 
from the CMB alone (Jaffe et al. 2006a,c). A Bayesian analysis of 
Bianchi Vll h models was performed by Bridges et al. (2007) using 
WMAP ILC data to explore the joint cosmological and Bianchi 
parameter space via Markov chain Monte Carlo sampling. The 
Q m -Q A degeneracy of Bianchi Vll h models was studied more thor- 
oughly and, although a similar best-fit Bianchi template was found 
in WMAP data, it was again determined that the parameters of the 
resulting Bianchi cosmology were inconsistent with standard con- 
straints. In a following study by Bridges et al. (2008) it was sug- 
gested that the CMB 'cold spot' (Vielva et al. 2004; Cruz et al. 
2006; Vielva 2010) could be driving the best-fit Bianchi compo- 
nent found in WMAP data. 

In this article we perform a definitive study of Bianchi Vllh 
cosmologies with WMAP data. We develop a Bayesian analysis 
technique capable of studying Bianchi models using partial- sky 
CMB observations, allowing us for the first time to study individual 
WMAP bands rather than the ILC map, which was not originally in- 
tended for cosmological analysis (Bennett et al. 2003). We sample 
the complete set of parameters describing the ACDM cosmology 
and the Bianchi Vllh model simultaneously, and perform the first 
rigorous study of the physically motivated scenario where the pa- 
rameters of the Bianchi model are coupled to the standard cosmol- 
ogy. To make comparisons with previous work, we also consider 
the non-physical scenario where the Bianchi parameters are decou- 
pled from the standard cosmological parameters. We employ the 
latest (and final) 9-year release of WMAP observations (previous 
studies have considered WMAP 1- and 3-year data only), and per- 
form a Bayesian model- selection analysis of Bianchi Vllh models, 




Figure 1. Simulated deterministic CMB temperature contributions in 
Bianchi Vllh cosmologies for varying x and £2 t otai (left-to-right £2 t otai e 
{0.10,0.30,0.95}; top-to-bottom x e {0.1,0.3,0.7, 1.5,6.0}). In these maps 
the spiral pattern typical of Bianchi Vllh induced temperature fluctuations 
is rotated from the South pole to the Galactic centre for illustrational pur- 
poses. The parameter values k = +1 and (a,/3, y) = (0°, -90°, 0°) were also 
set when generating these simulations. 

using nested sampling methods (Skilling 2004; Feroz & Hobson 
2008; Feroz et al. 2009), to determine whether WMAP data sug- 
gest we inhabit an anisotropic Bianchi Vll h universe instead of a 
standard isotropic ACDM universe. 

The remainder of this article is structured as follows. In 
Sec. 2 we briefly review Bianchi Vllh cosmologies. We describe 
the Bayesian analysis techniques that we develop to study full- and 
partial-sky CMB observations in Sec. 3. Although focus is given 
to Bianchi models, these techniques are generic and may be used 
to study other general anisotropic cosmologies. In Sec. 4 we ap- 
ply our analysis to WMAP 9-year data to constrain Bianchi Vll h 
cosmologies. Concluding remarks are made in Sec. 5. 



2 BIANCHI \ll h COSMOLOGIES 

In this section we review Bianchi Vll h cosmologies, focusing on 
their description; we defer technical details to Barrow et al. (1985) 
and Jaffe et al. (2006c). Note that we adopt the solutions to the 
temperature fluctuations induced in Bianchi Vll h models when in- 
corporating dark energy, as derived by Lasenby and also adopted 
by Bridges et al. (2007). 

Bianchi VII^ models contain a free parameter, usually denoted 
x, as first identified by Collins & Hawking (1973). This parameter 
is related to the h parameter of type Vll h models by 

where the total energy density Q tota i = n m + Q A . Physically, x is re- 
lated to the characteristic wavelength over which the principle axes 
of shear and rotation change orientation. Consequently, x acts to 
change the 'tightness' of the spiral-type CMB temperature contri- 
butions that are typical of Bianchi Vll h cosmologies. Bianchi Vll h 
models are also described by their shear (o-/H) , vorticity (q)/H)q, 
and handedness /c, where H is the Hubble parameter. The handed- 
ness parameter takes the values k = +1 and k - -1 for maps with 
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right- and left-handed spirals respectively. Vorticity is related to the 
other parameters by (Barrow et al. 1985) 

\hI 3 V2^n total \Hk' 

where equality of shear modes is assumed (cf. Jaffe et al. 2005). The 
amplitude of deterministic Bianchi VII /7 temperature maps may be 
characterised by either (<r/H)o or (q)/H)q since these parameters 
influence the amplitude of the map only and not its morphology. 
Since the CMB temperature fluctuations induced in Bianchi VII n 
models are anisotropic on the sky, the orientation of the CMB con- 
tribution may vary also, introducing three additional degrees-of- 
freedom. The orientation of a map representing the Bianchi VII/j 
CMB fluctuations is described by the Euler angles 1 (a,p, y), where 
for (a,fi,y) = (0°, 0°, 0°) the spiral pattern typical of Bianchi \ll h 
temperature fluctuations is centred on the South pole. To sum- 
marise, Bianchi Vll h models may be described by the parameter 
vector B = (n m , £2 A , x, (a>/H)o, a,j3,y) (note that we do not in- 
clude the handedness parameter k in B since this is used to distin- 
guish between left- and right-handed models). 

In the analysis performed herein the BIANCHI 2 2 code is used 
to simulate the temperature fluctuations induced in Bianchi VII/, 
models. Bianchi VII/ 7 models induce only large scale temperature 
fluctuations in the CMB and, consequently, the resulting Bianchi 
maps have a low band-limit both globally and azimuthally, i.e. 
in both i and m in spherical harmonic space; indeed, only those 
harmonic coefficients with m = +1 are non-zero (Barrow et al. 
1985). In the BIANCHI2 code Bianchi VII/, temperature fluctua- 
tions may be computed directly in either real or harmonic space. 
In the analysis performed herein we compute temperature fluctu- 
ations directly in harmonic space to avoid any pixelisation arte- 
facts and since the resulting temperature fluctuations can be rotated 
accurately and efficiently in harmonic space (see McEwen et al. 
2006) due to their low azimuthal band-limit. Examples of simu- 
lated Bianchi Yll h temperature maps are illustrated in Fig. 1 for a 
range of parameters. 



3 BAYESIAN ANALYSIS OF ANISOTROPIC 
COSMOLOGIES 

We describe in this section a generic methodology to perform a 
Bayesian analysis of anisotropic cosmologies using CMB observa- 
tions. Since we are motivated by the analysis of Bianchi VII/, cos- 
mologies, we consider a deterministic global template and pose the 
analysis in harmonic space (where Bianchi VII/, contributions can 
be computed accurately and rotated efficiently). The extension to 
alternative settings, such as non-trivial topologies (Niarchou et al. 
2003; Cornish et al. 2004; Dineen et al. 2005; Kunz et al. 2006; 
Bielewicz & Riazuelo 2009), is also highlighted. Firstly, we de- 
scribe a generic Bayesian analysis of anisotropic cosmologies, be- 
fore restricting this to the specific settings of full- and partial-sky 
CMB observations. Finally, we perform simulations to validate the 
implementation of the methodology presented. 



The active zyz Euler convention is adopted, corresponding to the rotation 
of a physical body in a fixed coordinate system about the z, y and z axes by 
y, and a respectively. 
2 http://www.jasonmcewen.org/ 



3.1 Bayesian analysis 

By B ayes' theorem we may write the posterior distribution of the 
parameters of our model of interest M, given data d, as 

P(e|.^ = p(rf|0 p Xf |M) ocP(.|e,M)P ( e|M), (3) 

P(d | M) 

where P(d | 0, M) is the likelihood, P(0 1 M) is the prior distribu- 
tion of the parameters of the model, and P(d | M) is the Bayesian 
evidence, which normalises the posterior distribution. The (unnor- 
malised) posterior distribution encodes our inferred knowledge of 
the parameters of the model, given the observational data. 

We consider both models that include a Bianchi VII/, contribu- 
tion in addition to a stochastic CMB component and those that do 
not. We consider open and flat cosmologies where the cosmologi- 
cal parameters are given by C = (A s , n s , r, Qbh 2 , £l c h 2 , Qa, £lk), 
where A s is the amplitude of the primordial power spectrum, n s 
is the scalar spectral index, r is the optical depth of reionisation, 
Q h h 2 is the physical baryon density, Q c h 2 is the physical cold dark 
matter density, Q A is the dark energy density, £l k is the curvature 
density and H - 100 h; for the flat case Q k = and we recover the 
standard six parameter model. For the models including a Bianchi 
VII/, contribution, we consider the physical case where the Bianchi 
parameters B and cosmological parameters C are coupled (i.e. 
where the cosmological density parameters shared by B and C 
are set equal) and also the non-physical case where they are not (the 
latter case is considered to make comparisons with previous work). 
The Bianchi Vll h parameters B are described in Sec. 2. 

For the case where a Bianchi VII/, template is embedded in 
a stochastic CMB background described by its power spectrum 
Q(0 C ), the likelihood is given by 

P(<* |©b,0c) oc 1 exp[-^ 2 (0c,0 B )/2] , (4) 
V|X(0c)| 

where 

X 2 (®c, ®b) = [d- b(® B )f X _1 (0 C ) [d - b(Q B )] , (5) 

and where #(0 B ) is the deterministic Bianchi VII/, template and 
X(0 C ) is the covariance matrix of the stochastic CMB component 
(and any noise component included in the model). Eqn. (4) and 
Eqn. (5) are written in a generic manner, where the likelihood may 
be given either in real or harmonic space. Since Bianchi VII/, tem- 
plates can be computed accurately and rotated efficiently in har- 
monic space, we specialise to a harmonic space representation, 
where the data and Bianchi template are given by their spherical 
harmonic coefficients: d = {d im } and Z>(0 B ) - {^(©b)}, respec- 
tively, considered up to the harmonic band-limit £ max . The covari- 
ance matrix X(0 C ) is then also specified in harmonic space but dif- 
fers depending on whether the full- or partial- sky setting is consid- 
ered (we consider each setting in turn in subsequent subsections). 
For the full- sky setting the covariance matrix is diagonal, whereas 
this is not the case for the partial-sky setting (assuming isotropic 
noise in both settings). 

We have so far considered the setting where a subdominant, 
deterministic template is embedded in a stochastic CMB back- 
ground. However, it is straightforward to extend the analysis to the 
setting where an anisotropic cosmology does not induce an embed- 
ded deterministic template, but rather alters the covariance structure 
of the stochastic CMB component; this is the case for non-trivial 
topologies (e.g. Niarchou et al. 2003; Cornish et al. 2004; Dineen 
et al. 2005; Kunz et al. 2006; Bielewicz & Riazuelo 2009). This set- 
ting can be handled in the current framework simply by setting the 
template component to zero and by replacing the covariance matrix 
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X(0 C ) with the appropriate covariance matrix for the anisotropic 
cosmology. The coordinate orientation is again arbitrary in this set- 
ting, hence the covariance matrix, or similarly the data, must be 
varied over all three-dimensional orientations as specified by the 
Euler angles. In this setting the^ 2 of the likelihood is written 



X 2 (®a) = et(a,p,y)X-\® K )d(a,p,y) , 



(6) 



where A is the full set of parameters of the anisotropic cosmol- 
ogy and d(a,/3, y) denotes the data rotated by Euler angles (a,J3, y). 
Since we focus on Bianchi cosmologies here, we consider only the 
likelihood with^f 2 given by Eqn. (5) in the remainder of this article. 

To determine whether the inclusion of a Bianchi VII/ 7 compo- 
nent better describes the data the Bayesian evidence is examined, 
given by 



E = P(d\M) = J d0P(rf|0,M)P(0|M). 



(7) 



Using the Bayesian evidence to distinguish between models natu- 
rally incorporates Occam's razor, trading off model simplicity and 
accuracy. The odds ratio between two models Mi and M 2 may be 
written in terms of the Bayesian evidences of the models (E\ and 
E 2 respectively) by 



P(M 2 \d) 



Ex P(MQ 
E 2 P(M 2 ) 



(8) 



In the absence of any prior information about a preferred model, 
i.e. when P(Mi) = P(M 2 ), the Bayes factor given by the ratio of 
Bayesian evidences Ei/E 2 is thus identical to the ratio of the model 
probabilities given the data. Consequently, the Bayes factor may 
used to distinguish models in this setting. 

The Jeffreys scale (Jeffreys 1961) is often used as a rule-of- 
thumb when comparing models via their Bayes factor. The log- 
Bayes factor A\nE = \n(Ei/E 2 ) represents the degree by which 
model Mi is favoured over model M 2 , assuming the models are 
equally likely a priori. On the Jeffreys scale log-Bayes factors are 
given the following interpretation: < A\nE < 1 is regarded as 
inconclusive; 1 < A\nE < 2.5 as significant; 2.5 < A\nE < 5 as 
strong; and A\nE ^ 5 as conclusive (without loss of generality we 
have assumed E\ ^ E 2 ). For reference, a log-Bayes factor of 2.5 
corresponds to odds of approximately 1 in 12, while a factor of 5 
corresponds to odds of approximately 1 in 150. 



3.2 Full-sky analysis 

In the full- sky setting the covariance matrix is given by 
X(0 C ) = C(0 C ), where C(0 C ) is the diagonal CMB covariance ma- 
trix with entries given by the power spectrum Q(0 C ) on its diag- 
onal. A Bayesian analysis in this setting was first considered by 
Bridges et al. (2007). To expose the detail of the analysis we make 
the likelihood explicit in this setting. The likelihood, in terms of 
the spherical harmonic coefficients of the data d im and Bianchi Vll h 
template b^ n (S B )^ is given by 



P({d (m ) | 0b, C ) oc PI - exp - 
l e J VQ(®c) I 



[dm - MOn)r 
Q(0 C ) 



A 2 / \d tm -b {m (& B f \ 

1 1 ^7zr\ ex P1 — \ ■ (9) 



Q(©c) 



Q(©c) 



In practice, it is more convenient numerically to work with the log- 
likelihood function, given by 



|0B,0C)]«2| ( 

(=0 v 



ln[P({fiU | ©b, 0c)] { (2£ + l)ln[Q(0 c )] 

t 

[d„ - b eo (® B )] 2 2 



Q(©c) 



Q(0 C ) 



J]\d em - b em (®B)\ 2 ) ■ (10) 

m=1 / 



In the case of a zero Bianchi component b(® B ) ~ 0> Eqn. (10) 
reduces to the log-likelihood function used commonly to compute 
parameter estimates from the power spectrum estimated from CMB 
data (e.g. Verde et al. 2003). 



3.3 Partial-sky analysis 

In the partial- sky setting a mask is applied in real space to excise 
contaminated regions of the data (due to point- source or Galactic 
contamination, for example). See, for example, the mask applied to 
simulated CMB data in Fig. 4(a). To account in harmonic space for 
a mask applied in real space, Gaussian masking noise m is added to 
the data d in real space: 



s(oj) = d{(jj) + m((jj) , 



en) 



where cj denotes angular coordinates on the sphere. 

The masking noise is chosen to be zero-mean and large in the 
masked region of the data (i.e. zeros of the mask), and zero else- 
where. Consequently, the masking noise is anisotropic over the sky 
but may be chosen to be uncorrelated, and may thus be defined by 
its covariance 



(m(cx)i)rn(cx)j)) = 6^0^(0)1) , 



(12) 



where is Kronecker delta symbol, aji denotes the angular coor- 
dinate of pixel i, and the variance of the noise for pixel i is given 
by 



o- m (0Ji) ■ 



0, 



oji e } 



(13) 



where E 2 is a constant masking noise variance. We adopt the 
HEALPix (Gorski et al. 2005) equal-area pixelisation of the sphere, 
where the area of each pixel is Q pix , in order to avoid dealing with 
pixels of differing areas. Here M denotes the masked region of the 
sky S 2 and S 2 \M denotes the remaining region. By synthetically 
adding masking noise that is much larger than the original data in 
the masked region of the sky, we effectively marginalise the pixel 
values of the data in this region. Consequently, only the pixel val- 
ues of the data outside of the masked region have a large influence 
on the final analysis. We refer to the a particular masking noise 
realisation as a noisy mask. An example noisy mask is shown in 
Fig. 4(d). 

The noisy mask must be accounted for in the analysis. The ef- 
fective signal s to be analysed now includes the original data d and 
the noisy mask m; thus its covariance is the sum of the covariance 
of the original data and the noisy mask. The covariance of the noisy 
mask is given in harmonic space by 



M, 



: {m €m m\, m , > ^ cr 2 (aJi)Y* em (aJi) Y Vm > (aj t ) Q pix 2 , (14) 



where we have applied Eqn. (12) and again consider signals and 
their covariance structure in harmonic space (since the analysis of 
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Bianchi Vll h cosmologies is most accurately and efficiently per- 
formed in harmonic space). In Eqn. (14) we approximate the spher- 
ical harmonic coefficients of the mask using discrete quadrature: 

m tm = I 6Q(cj) m(oj) Y* m (aj) V Q pix m{oj { ) Y* m (oji) ; (15) 

note that this is necessarily an approximation since the noisy 
mask is not band-limited and we adopt the HEALPix (Gorski 
et al. 2005) sampling scheme, which does not admit exact quadra- 
ture. The co variance of the effective data s is thus given by 
X(0 C ) = C(0 C ) + M, where C(0 C ) is again the diagonal CMB co- 
variance matrix and M is the non-diagonal noisy mask covariance 
matrix given by Eqn. (14). 

To summarise, the partial-sky analysis proceeds as follows. 
Firstly a zero-mean, Gaussian noisy mask realisation is con- 
structed via Eqn. (12) and Eqn. (13), and its covariance struc- 
ture is estimated by Eqn. (14). The effective signal under analy- 
sis s is constructed by summing the original data and the noisy 
mask realisation by Eqn. (11). The Bayesian analysis described 
in Sec. 3.1 is performed, where the data under analysis is re- 
placed with the effective signal s, with covariance matrix given by 
X(0 c ) = C(0 c ) + M. 3 

3.4 Implementation and validation 

We have developed the ANIC0SM0 4 code to perform a Bayesian 
analysis of anisotropic cosmologies. The Bayesian analysis de- 
scribed in this section is implemented in ANIC0SM0, both in the full- 
and partial-sky settings. We adopt the MultiNest 5 code (Feroz & 
Hobson 2008; Feroz et al. 2009) to explore the posterior distribu- 
tion and to compute Bayesian evidence values, both via nested sam- 
pling (Skilling 2004; Feroz & Hobson 2008; Feroz et al. 2009). The 
CAMB 6 code (Lewis et al. 2000) is used to compute the CMB power 
spectrum Q(0 C ) for various cosmological parameters C . 

At present the ANIC0SM0 code is specialised to the study 
of Bianchi VII^ cosmologies, where the likelihood is computed 
through Eqn. (4) and Eqn. (5). However, ANIC0SM0 may be triv- 
ially extended to handle the x 2 defined by Eqn. (6) and to study 
other anisotropic cosmologies, such as non-trivial topologies; this 
will be the focus of future work. 

In the remainder of this subsection we perform simulations to 
validate the Bayesian analysis method described here and its im- 
plementation in the ANIC0SM0 code. In both the full- and partial- 
sky settings, we simulate a CMB map and synthetically embed 
a simulated Bianchi VII^ component. To approximate the finite- 
resolution, noisy measurements of the CMB made by the WMAP 
experiment (specifically those made in the highest-resolution, 94 
GHz W band; Bennett et al. 2012), a Gaussian beam with full- 
width-half-maximum (FWHM) 13.2' and isotropic noise with 
power spectrum Nt = 0.02 {jiK) 2 are applied to the simulations. 

3 In practice the masked region of the data to be excised will include un- 
known contamination, rather than a CMB contribution only. Thus, the orig- 
inal data d should be masked before the procedure outlined here is per- 
formed. In this setting the covariance matrix X(@c) = C(©c) + M is neces- 
sarily an approximation. However, provided the noise variance is chosen 
to be considerably larger than the expected CMB contribution, the approxi- 
mation will be very accurate. 

4 The AN I COSMO code will soon be made publicly available from 
http : //www . j asonmcewen . org/. 

5 http : //www . mrao . cam . ac . uk/sof tware/multinest/ 

6 http://camb.info/ 
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Figure 2. Full-sky simulation at £ max = 512. A Bianchi VII/j component 
with a large amplitude is simulated (panel b) and embedded in a standard 
stochastic CMB map (panel a), yielding a resultant CMB plus Bianchi map 
(panel c) to which a beam and noise are applied. The full-sky Bayesian anal- 
ysis described in Sec. 3.2 is applied to recover a MAP estimated Bianchi 
component (panel d). The recovered Bianchi component (panel d) accu- 
rately estimates the embedded component (panel b). 

To account for the beam and noise in the Bayesian analysis we 
simply map Q(@ c ) — > b 2 € Q(0 C ) + Nt when computing the CMB 
covariance matrix and apply the beam to the Bianchi Nll h template 
by mapping bt m (S B ) — > bt ^(©bX where b 2 and Nt are the beam 
and noise power spectra, respectively. In these simulations we con- 
sider the model where the Bianchi parameters B and cosmologi- 
cal parameters C are coupled. Posterior distributions and evidence 
values are then recovered by ANIC0SM0 for models including and 
excluding a coupled Bianchi component in addition to the stan- 
dard stochastic CMB component. For the purpose of validating the 
implementation, we consider a Bianchi template with a relatively 
large amplitude of (oj/H) = 2 x 10~ 9 . 
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Figure 3. Posterior distributions corresponding to the full-sky simulation at £ mSLX = 512 shown in Fig. 2. The underlying parameter values of the simulation 
are indicated by dashed vertical lines. Both cosmological and Bianchi VII/ 7 parameters are estimated accurately, except for the optical depth of reionisation, r, 
which is not unexpected when using CMB temperature data alone. For model comparison purposes an unconstrained r does not pose any great concern. 



Full- sky simulations are shown in Fig. 2. Since the covari- 
ance matrix X(0 C ) is diagonal in this setting, the analysis may 
be run at high-resolution; hence, we set £ mSLX = 512 in this simu- 
lation. The maximum-a-posteriori (MAP) estimated Bianchi VII/j 
component recovered from the peak of the posterior distribution 
P(0 1 d, M) is shown in Fig. 2(d) and is clearly an accurate estimate 
of the embedded component shown in Fig. 2(b). The marginalised 
posterior distributions recovered for each parameter are shown in 
Fig. 3. Both cosmological and Bianchi Nll h parameters are es- 
timated accurately, except for the optical depth of reionisation r 
which is not unexpected when using CMB temperature data alone, 
r is also unconstrained when considering the model absent of a 
Bianchi component, thus for model comparison purposes an un- 
constrained r does not pose any great concern. Indeed, for this 
simulation the model including a Bianchi component is favoured 
conclusively with a log-Bayes factor of Aln£ ~ 50. 



Partial-sky simulations are shown in Fig. 4. Since the covari- 
ance matrix X(0 C ) is non-diagonal in this setting, the analysis is 
considerably more computationally demanding than the full- sky 
setting; hence, we set £ max = 32 in this simulation. The MAP esti- 
mated Bianchi YIl h component is shown in Fig. 4(e) and clearly 
is an accurate estimate of the embedded component shown in 
Fig. 4(b) (note that the embedded Bianchi component was rotated 
such that the majority of its structure does not lie in the masked 
region). A noisy mask variance of = 100 (mK) 2 was adopted 
since numerical tests showed this effectively marginalised masked 
pixels in the data, while increasing S 2 further had little impact. The 
marginalised posterior distributions recovered for each parameter 
are shown in Fig. 5. Both cosmological and Bianchi VII^ parame- 
ters are estimated reasonably accurately, however since less data is 
now available due to the lower £ max and masking, the marginalised 
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Figure 5. Posterior distributions corresponding to the partial- sky simulation at ^ max = 32 shown in Fig. 4. The underlying parameter values of the simulation 
are indicated by dashed vertical lines. Both cosmological and Bianchi VII^ parameters are estimated reasonably well, however since less data is now available 
due to the lower ^ max and masking, the marginalised posterior distributions are not as accurate or well constrained as the full-sky setting, as is to be expected. 
For model comparison purposes this does not pose any great concern. 



posterior distributions are not as well constrained as the full- sky 
setting, as is to be expected. 

Reducing the maximum multipole considered has a large im- 
pact on the ACDM parameters determined by the acoustic peak 
positions and heights, but, thanks to the effective low-^ band-limit 
of the Bianchi Vll h signal, affects the Bianchi Yll h parameters 
much more mildly. This is clearly evident from the comparative 
behaviour of the posteriors for (a)/H) and, e.g., Q c h 2 or Q h h 2 . For 
model comparison purposes, this does not therefore pose any great 
concern. Indeed, for this simulation the model including a Bianchi 
component is also favoured conclusively with a log-Bayes factor 
of Aln£ ~ 50. Nevertheless, by reducing £ mSLX we clearly discard 
a great deal of cosmologically interesting information. This can be 
recovered by calculating a high-^ likelihood assuming the Bianchi 
VII/j contribution is zero for £ > ^ max , using a conjugate-gradient- 



based implementation (Smith et al. 2007; Bennett et al. 2012) of 
the optimal, unbiased power- spectrum estimator (Tegmark 1997) 
in the presence of a mask. The two likelihoods can then be multi- 
plied, assuming that the high- and low-^ data are independent. The 
implementation of such an algorithm for the simulations used here 
is unnecessary, as they are intended solely to validate the ANI COSMO 
code. When considering the WMAP data, we will use the purpose- 
built WMAP likelihood code 7 (Bennett et al. 2012), which imple- 
ments the scheme described above, to calculate the high-^ likeli- 
hood and hence maximise the impact of the WMAP data. 

The rationale for the simulations performed here is to vali- 



7 Available for download from http://lambda.gsfc.nasa.gov/ 
product/map/dr5/likelihood_get . cfm. 
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Figure 4. Partial-sky simulation at £ max = 32. A Bianchi VII^ component 
with a large amplitude is simulated (panel b) and embedded in a standard 
stochastic CMB map (panel a), yielding a resultant CMB plus Bianchi map 
(panel c) to which a mask, beam and noise are applied (the Bianchi com- 
ponent is rotated such that the majority of its structure does not lie in the 
masked region). A noisy mask (panel d) is added to the CMB plus Bianchi 
map and the partial-sky Bayesian analysis described in Sec. 3.3 is applied 
to recover a MAP estimated Bianchi component (panel e). The recovered 
Bianchi component (panel e) accurately estimates the embedded compo- 
nent (panel b). 



date the implementation of the ANIC0SM0 code, which has been 
demonstrated effectively in both the full- and partial- sky settings. 
Furthermore, although the details are not presented here, additional 
simulations were performed to validate the partial-sky setting. Two 
CMB maps were simulated with different cosmologies. A single 
hybrid map was then constructed with the Northern hemisphere 
given by the CMB simulated from the first cosmology and with 
the Southern hemisphere given by the CMB simulated from the 
second cosmology. ANIC0SM0 was applied to recover cosmologi- 
cal parameters from this hybrid map, once masking the Northern 
hemisphere and once masking the Southern hemisphere. In each 
of these tests the correct cosmology of the unmasked hemisphere 
was recovered, demonstrating that the partial- sky Bayesian analy- 
sis outlined in Sec. 3.3 effectively marginalises the masked pixels 
of the data. This test and the simulations described previously pro- 
vide a strong validation of the ANIC0SM0 code in both the full- and 
partial-sky settings. 



4 ANALYSIS OF WMAP OBSERVATIONS 

We analyse WMAP 9-year data for evidence of Bianchi YIl h cos- 
mologies, performing both the full- and partial-sky Bayesian anal- 
yses described in the preceding section. Firstly, we describe the 
specific full- and partial- sky WMAP data used and the cosmologi- 
cal models considered. We then present the results of our Bayesian 
analysis of Bianchi Yll h models and place robust constraints on the 
vorticity of the Universe. 



4.1 Data 

Previous searches for Bianchi VII/j components embedded in 
WMAP data (e.g. Bridges et al. 2007) used full-sky data in the 
form of the ILC map to estimate the parameters of the Bianchi 
VII/ 7 model. While this greatly simplifies the form of the covari- 
ance matrix employed in the likelihood, the ILC map was not rec- 
ommended for cosmological analysis (Eriksen et al. 2004) as it 
contains considerable foreground residuals, especially within the 
Galactic plane, and has a complex pixel noise structure. While these 
effects were ignored in previous studies (e.g. Bridges et al. 2007) 
since they are sub-dominant on the large scales of interest when 
studying Bianchi models, a better approach is to analyse partial- 
sky observations of individual WMAP frequency bands. 

Our formalism, presented in Sec. 3.3, allows partial-sky data 
to be analysed for the first time in a statistically rigorous manner, 
minimising the residual contamination in the data and ensuring our 
conclusions are robust. We analyse the foreground-reduced WMAP 
9-year W-band map (Bennett et al. 2012), masked with the conser- 
vative 9-year KQ75 sky cut (Bennett et al. 2012) to excise residual 
contamination. We select W-band data since this band, along with 
the V-band (which has lower resolution), suffers from the least fore- 
ground contamination (Bennett et al. 2012). The resulting data are 
shown in Fig. 6(a). In order to draw comparisons with previous 
studies we also analyse the full-sky WMAP 9-year ILC map (Ben- 
nett et al. 2012) shown in Fig. 6(b), under the proviso that any cos- 
mological conclusions drawn from the ILC analysis must be treated 
with care. 

Since Bianchi models have a low harmonic band-limit, we use 
only low-^ WMAP data in computing the Bianchi likelihood de- 
fined by Eqn. (4); specifically, we compute contributions to the 
likelihood of Eqn. (4) for i < ^ ax = 32. For the partial-sky set- 
ting, where a non-diagonal covariance matrix must be inverted for 
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each likelihood evaluation, a low band-limit ^ ax is also convenient 
to reduce the computational cost of the analysis. The WMAP data 
so far described (the ILC map or the KQ75-masked W-band map 
in the full- and partial- sky settings respectively) are used for the 
Bianchi likelihood evaluation. As mentioned previously, we aug- 
ment the low-^ Bianchi likelihood with the standard high-^ WMAP 
likelihood (Bennett et al. 2012) for £ > ^ ax = 32; this likelihood 
function makes use of all WMAP temperature observations, as well 
as accurate beam models. The final log-likelihood is thus the sum 
of the \ow-£ Bianchi and high-^ 9-year WMAP log-likelihood con- 
tributions (since the low- and high-^ data are essentially indepen- 
dent, the log-likelihood contributions can effectively by summed). 
As highlighted in Sec. 3.4, by including the high-^ WMAP likeli- 
hood we are able to constrain cosmological parameters to greater 
precision. 

We incorporate approximate WMAP noise and beam effects 
in the \ow-£ Bianchi likelihood. Since the Bianchi likelihood is 
specified in harmonic space (due to the efficiency and accuracy for 
which Bianchi VII/ 7 models can be handled in harmonic space), 
we approximate the anisotropic WMAP noise by isotropic noise 
(which may be handled easily in harmonic space due to its di- 
agonal covariance structure). We assume isotropic white noise 
specified by its power spectrum N{ = 0.015 (fiK) 2 (computed by 
Nt = ^ P i x cr^/median(A^ obs ), where N ohs is the W-band observation- 
count map and cr = 6544 fiK is the W-band pixel noise level) and 
a Gaussian beam with FWHM of 13.2' in order to approximate 
W-band observations accurately. As the individual WMAP bands 
are smoothed to a common resolution of 1° when creating the ILC 
(Bennett et al. 2012), we substitute a Gaussian beam of FWHM 
1° when considering ILC data. Since we restrict to low-^ data for 
the Bianchi likelihood {£ < ^ ax = 32), where WMAP noise and 
beam effects are sub-dominant, these approximations are accurate. 
For the high-^ contribution to the likelihood, the official WMAP 9- 
year likelihood is used, where noise and beams are modelled very 
accurately. 

4.2 Models 

We consider three scenarios in this study: one physically motivated, 
and two phenomenological but motivated by previous analyses. 

In the first (and only physically motivated) scenario, the 
Bianchi and cosmological parameters are coupled; i.e. the matter 
and dark-energy densities of the Bianchi and standard cosmologi- 
cal models are shared and thus identical. Since Bianchi VII /7 models 
are open (with flat models as the limiting case) we label this model 
the open-coupled-Bianchi model. When quoting final conclusions 
about Bianchi VII/j cosmologies we will refer to this model since it 
is a physically well-motivated and consistent model. 

It is also interesting to consider Bianchi CMB contributions 
arising from phenomenological models in which the Bianchi pa- 
rameters are decoupled from the standard cosmological parame- 
ters; indeed, it is necessary to do so in order to compare to existing 
results, as this is the only approach employed in all previous anal- 
yses. In this non-physical scenario we consider a standard ACDM 
cosmology and a Bianchi VII/ 7 component with decoupled param- 
eters {i.e. the total energy density of the standard cosmology and 
Bianchi cosmology may differ and are effectively separate param- 
eters). As the ACDM cosmologies considered in previous analyses 
{e.g. Jaffe et al. 2005; Bridges et al. 2007) were flat, we similarly re- 
strict our attention to flat ACDM components, and label this model 
the flat-decoupled-Bianchi model. Note that in this model we allow 
all ACDM parameters to vary. 




(a) Foreground-reduced W-band map with KQ75 mask applied 




600 fiK 



(b) ILC map 
Figure 6. WMAP 9-year data analysed. 



Finally, we also consider the situation where the majority of 
the cosmological parameters are fixed and are not fitted simultane- 
ously with the Bianchi Yll h parameters. Following Bridges et al. 
(2007), we do allow the amplitude of the primordial power spec- 
trum, A j, to vary in this model, so that the amplitude of standard 
stochastic CMB temperature fluctuations is allowed to vary as we 
fit for an embedded Bianchi component. The remaining cosmolog- 
ical parameters are fixed at their values as constrained by WMAP 
9-year observations, baryon acoustic oscillations and supernovae 
observations (Hinshaw et al. 2012). This is again a decoupled, non- 
physical scenario, and we label this model the fixed-decoupled- 
Bianchi model, where 'fixed' indicates that the standard cosmo- 
logical parameters are essentially fixed. 

For each of the three models discussed above, we consider 
models where a Bianchi VII/ 7 component is, and is not, included. 
Moreover, for the case where a Bianchi component is included, 
we consider both left-handed (LH) and right-handed (RH) Bianchi 
models, since the handedness of the coordinate system is also free 
in Bianchi VII/ 7 models. Consequently, we recover nine different 
models that we analyse. 

The prior ranges adopted on the ACDM and Bianchi parame- 
ters in this analysis are shown in Table 1 , and are chosen conserva- 
tively to reflect our weak prior knowledge of the Bianchi models. 
The priors are uniform in all parameters of interest apart from the 
prior on the power spectrum amplitude, A s , which is uniform in 
log A s , and the prior on the Euler angle J3, which is uniform in sin 
p. The priors on the Bianchi densities, and Q^, are applied only 
in the decoupled models; in the coupled model these parameters are 
set by the sampled ACDM densities according to = £l h + Q c 
and Ql = Q A . 
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Table 1. The prior ranges used for the ACDM and Bianchi VII/j parameters. 



Parameter 


Prior Range 


A s 


[1,5] xlO" 9 


n s 


[0.8,1.1] 


T 


[0.082,0.092] 


£l h h 2 


[0.005, 0.05] 


O //,2 


[U.UJ, U. jJ 




[0.5,0.9] 


Q k 


[0.001,0.2] 


X 


[0.01, 1] 


(o)/H) 


[0, 1] x 10" 10 


a 


[0, 360]° 


P 


[0,180]° 


7 


[0, 360]° 




[0,0.99] 




[0,0.99] 



Table 2. Log-Bayes factors computed for different models. The Bayes 
factor for each model is computed relative to the corresponding model 
that does not include a Bianchi component, where a positive Bayes factor 
favours the model that does include a Bianchi component. 



Model Full-sky ILC Masked W-band 



Open-coupled-Bianchi (LH) 


-0.3 ± 0.2 


0.0 


± 


0.2 


Open-coupled-Bianchi (RH) 


-0.3 ± 0.2 


0.1 


± 


0.2 


Flat-decoupled-Bianchi (LH) 


1.1 ±0.2 


0.1 


± 


0.2 


Flat-decoupled-Bianchi (RH) 


0.1 ± 0.2 


0.1 


± 


0.2 


Fixed-decoupled-Bianchi (LH) 


1.7 ±0.1 


0.6 


+ 


0.1 


Fixed-decoupled-Bianchi (RH) 


0.6 ±0.1 


0.3 


± 


0.1 



4.3 Results 

We use WMAP 9-year data, as specified in Sec. 4.1, to study the 
Bianchi Yll h models described in Sec. 4.2, using both the full- and 
partial-sky Bayesian analysis techniques presented in Sec. 3. We 
first consider the Bayes factors computed for different models, be- 
fore studying parameter estimates, best-fit maps, and finally vortic- 
ity bounds. 



4.3.1 Bayesian evidence 

The log-Bayes factors computed for the different models consid- 
ered are shown in Table 2. Neither left- nor right-handed open- 
coupled-Bianchi models are favoured by WMAP 9-year data, for 
either full- sky ILC data or partial- sky W-band data. Since this is 
the only physical model studied, we may already conclude that 
WMAP data do not favour a Bianchi Yll h cosmology over the stan- 
dard ACDM cosmology. 

Nevertheless, we consider the other phenomenological mod- 
els discussed in Sec. 4.2 where the Bianchi parameters are de- 
coupled from the standard cosmological parameters, since these 
are the models that have been studied in previous analyses. For 
both the left-handed, flat-decoupled-Bianchi and fixed-decoupled- 
Bianchi models, we find evidence in favour of a Bianchi compo- 
nent (classified as significant on the Jeffreys scale), when analysing 
full- sky ILC data. However, when analysing masked W-band data 
we find no evidence for these models, suggesting that the evidence 
found in the full-sky setting is driven predominantly by ILC data 
near the Galactic plane. No evidence is found for the correspond- 
ing right-hand models. 



4.3.2 Parameter estimates 

Although only left-handed, decoupled-Bianchi models show evi- 
dence for a Bianchi component when analysing full- sky ILC data, 
for completeness we show posterior distributions of the parameters 
of all left-handed models, using both full-sky ILC data and masked 
W-band data. We focus on left-handed models since neither dataset 
shows evidence for any right-handed models. 

The posterior distributions of the parameters of the open- 
coupled-Bianchi model are shown in Fig. 7 for full- and partial- 
sky data. The standard cosmological parameters are relatively well 
constrained thanks to the high-^ WMAP likelihood (except for the 
optical depth to reionisation, r, as expected as only CMB temper- 
ature data are used). Since the constraints on the standard cosmo- 
logical parameters are driven largely by the high-^ WMAP likeli- 
hood and are similar for all models, we do not comment on these 
further when discussing subsequent models. The Bianchi parame- 
ters, however, are typically poorly constrained for this model, es- 
pecially for partial-sky W-band data. Furthermore, the Bianchi vor- 
ticity (o)/H)o, which traces the amplitude of the Bianchi-induced 
CMB temperature fluctuations, is peaked near zero. These poste- 
rior plots agree with the conclusion from the Bayesian evidence 
that WMAP 9-year data do not favour Bianchi Yll h cosmologies. 

The posterior distributions of the parameters of the flat- 
decoupled-Bianchi model are shown in Fig. 8 for full- and partial- 
sky data. For the full-sky analysis, the Bianchi parameters are rela- 
tively well-constrained (except for due to the Qj^-Q^ degener- 
acy of Bianchi VII/ 7 models), agreeing with the previous inference 
from the Bayesian evidence that a Bianchi component is favoured 
in this setting. Bianchi parameters are not well constrained for 
partial-sky data and, moreover, the Bianchi vorticity (a>/H)o is 
peaked near zero, agreeing with the conclusion that partial-sky data 
do not favour a Bianchi component. 

The posterior distributions of the parameters of the fixed- 
decoupled-Bianchi model are shown in Fig. 9 for full- and partial- 
sky data. The conclusions drawn from the posterior distributions 
of the parameters of this model are identical to the flat-decoupled- 
Bianchi model: essentially, the Bianchi parameters are well con- 
strained when full- sky ILC data are used, but not when partial- sky 
data are used. Again, these findings agree with conclusions drawn 
previously using the Bayesian evidence. 

Finally, we summarise the parameter estimates made from the 
joint posterior distribution P(0 | d, M) for each model in Table 3 
and Table 4, for partial- and full-sky data respectively. We compute 
both maximum a-posteriori (MAP) parameter estimates from the 
peak of the joint posterior distribution and estimates from the mean 
of the marginalised posterior distribution for each parameter. One- 
standard-deviation errors are also quoted for the mean-posterior pa- 
rameter estimates. 



4.3.3 Best-fit maps 

For the two models where the Bayes factor provides significant evi- 
dence in support of a Bianchi component (see Table 2), namely the 
non-physical flat-decoupled-Bianchi and fixed-decoupled-Bianchi 
models, we plot the best-fit embedded Bianchi maps in Fig. 10(a) 
and Fig. 10(b), respectively. Best- fit maps are computed from the 
MAP parameter estimates contained in Table 3. These best-fit maps 
are very similar to those found in previous releases of WMAP data 
(Jaffe et al. 2005; Bridges et al. 2007). 
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Figure 7. Posterior distributions for parameters of the left-handed open-coupled-Bianchi model computed from full-sky ILC data (solid blue curves) and 
KQ75-masked W-band data (dashed red curves). 



4.3.4 Vorticity bounds 

We conclude our analysis of WMAP 9 -year data by reiterating that 
only the open-coupled-Bianchi model is physical. The other mod- 
els considered are non-physical since the Bianchi parameters are 
decoupled from the standard cosmological parameters: the Bianchi 
component can therefore require parameters, for example, a total 
energy density, that are incompatible with the stochastic CMB com- 
ponent. For the physical open-coupled-Bianchi model no evidence 
for a Bianchi component is found in either full-sky ILC or partial- 
sky W-band data. In this physical model we can therefore place 
constraints on the vorticity of Bianchi VII/ 7 cosmologies. For full- 
sky ILC data we recover the constraint (a)/H) < 8. 1 x 10~ 10 at 95% 
confidence, for both left- and right-handed models. For partial- sky 
W-band data we recover the constraint (qj/H) < 8.1 x 10~ 10 for the 
left-handed model and the constraint (a)/H) < 8.6 x 10~ 10 for the 
right-handed model, both at 95% confidence. These constraints on 



the global anisotropy of our Universe are placed in the context of a 
well-motivated, physical model, namely Bianchi VIIh cosmologies, 
using robust Bayesian statistical methods. 



5 CONCLUSIONS 

We have performed a definitive analysis of Bianchi VIIh cosmolo- 
gies with WMAP 9-year temperature data. We have studied full-sky 
ILC data and masked W-band data, using Bayesian analysis tech- 
niques developed for the full- and partial-sky settings. Three dif- 
ferent models were studied, namely the physically motivated open- 
coupled-Bianchi model, where the Bianchi and standard cosmo- 
logical parameters are coupled, and the flat-decoupled-Bianchi and 
fixed-decoupled-Bianchi models, which are non-physical since the 
Bianchi and standard cosmological parameters are decoupled. Al- 
though we have focused here on Bianchi models, our analysis tech- 
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Figure 8. Posterior distributions for parameters of the left-handed flat-decoupled-Bianchi model computed from full-sky ILC data (solid blue curves) and 
KQ75-masked W-band data (dashed red curves). 



niques may be easily extended to study other anisotropic cosmolo- 
gies such as non-trivial topologies; such models will be the focus 
of future research. 

For the non-physical decoupled models, we find Bayesian ev- 
idence favouring the inclusion of a left-handed Bianchi component 
when analysing full-sky ILC data. The resulting best-fit Bianchi 
maps found in WMAP 9-year data are similar to those found in 
previous releases of WMAP data (Jaffe et al. 2005; Bridges et al. 
2007). However, when studying these models with masked W-band 
data we find no evidence for a Bianchi component, suggesting that 
data near the Galactic plane provide a large contribution to the evi- 
dence found in the full- sky setting. 

For the physical Bianchi Vllh model, we find no evidence 
for the inclusion of a Bianchi component in either full- or partial- 
sky WMAP data. Since this is a well-motivated, physical model 
we can state definitively that WMAP data do not favour Bianchi 



Vllh cosmologies over ACDM. However, neither is it possible 
to conclusively discount Bianchi Vll h cosmologies in favour of 
ACDM cosmologies. We constrain the vorticity of Bianchi Vllh 
cosmologies, and hence the large-scale anisotropy of the universe, 
at (o)/H) < 8.6 x 10" 10 with 95% confidence. 

Bianchi Vll h cosmologies induce a global anisotropic CMB 
temperature contribution that is inherently low-resolution, so the 
higher resolution of forthcoming Planck (Planck collaboration 
2005) data as compared to WMAP is not in itself relevant to en- 
hanced constraints on these models. However, Planck's ability to 
remove large-scale astrophysical foregrounds will result in CMB 
observations of unprecedented precision over almost the entire sky. 
An analysis of Bianchi Vll h cosmologies using Planck temperature 
data will thus be very informative. Moreover, polarised Planck ob- 
servations will also have the potential to further constrain Bianchi 
Yll h cosmologies. 
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Figure 9. Posterior distributions for parameters of the left-handed fixed-decoupled-Bianchi model computed from full-sky ILC data (solid blue curves) and 
KQ75-masked W-band data (dashed red curves). 



Table 3. Parameter estimates recovered for various left-handed models from full-sky ILC data. Note that some of these models are not favoured by the Bayesian 
evidence and some parameters are not well constrained. We nevertheless show all parameter estimates for completeness. 



Parameter Open-coupled-Bianchi Flat-decoupled-Bianchi Fixed-decoupled-Bianchi 





MAP 


Mean 


MAP 


Mean 


MAP 


Mean 


A s 


2.21 x 10" 9 


(2.19 ± 0.03) x 10" 9 


2.17 x 10" 9 


(2.19 ± 0.04) x 10" 9 


2.219 x 10" 9 


(2.217 ± 0.009) x 10" 


n s 


0.98 


0.98 ± 0.01 


0.98 


0.98 ± 0.01 






T 


0.084 


0.087 ± 0.003 


0.084 


0.087 ± 0.003 






£l h h 2 


0.023 


0.0226 ± 0.0004 


0.023 


0.0225 ± 0.0004 






Q c h 2 


0.12 


0.112 ±0.004 


0.11 


0.112 ±0.004 






Q A 


0.72 


0.75 ± 0.03 


0.74 


0.73 ± 0.02 








0.002 


0.007 ± 0.005 










n£ 


0.27 


0.24 ± 0.03 


0.3 


0.3 ± 0.2 


0.4 


0.3 ±0.1 


^A 


0.72 


0.75 ± 0.03 


0.3 


0.3 ± 0.2 


0.2 


0.3 ± 0.2 


X 


0.3 


0.4 ± 0.3 


0.7 


0.6 ± 0.2 


0.6 


0.6 ± 0.2 


(o>/H)o 


8 x 10" 10 


(3 ± 2) x 10" 10 


10 x 10" 10 


(6 ± 3) x 10" 10 


9 x 10" 10 


(6 ± 2) x 10" 10 


a 


2.5° 


181.8° ± 100.1° 


41.3° 


88.8° ± 90.0° 


41.0° 


71.2° ±73.7° 




122.5° 


87.9° ± 35.7° 


27.2° 


50.8° ± 36.4° 


27.6° 


43.6° ± 30.9° 


J 


242.5° 


189.7° ± 100.9° 


319.8° 


259.9° ±90.1° 


314.1° 


277.7° ± 75.6° 



Table 4. Parameter estimates recovered for various left-handed models from KQ75-masked W-band data. Note that these models are not favoured by the 
Bayesian evidence and some parameters are not well constrained. We nevertheless show all parameter estimates for completeness. 



Parameter 


Open-coupled-Bianchi 


Flat-decoupled-Bianchi 


Fixed-decoupled-Bianchi 




MAP 


Mean 


MAP 


Mean 


MAP 


Mean 


A s 


2.17 x 10" 9 


(2.18 ± 0.04) x 10" 9 


2.19 x 10" 9 


(2.18 ± 0.04) x 10" 9 


2.226 x 10" 9 


(2.217 ± 0.008) x 10" 9 


n s 


0.98 


0.99 ± 0.01 


0.98 


0.98 ± 0.01 






T 


0.084 


0.087 ± 0.003 


0.084 


0.087 ± 0.003 






£l h h 2 


0.023 


0.0228 ± 0.0005 


0.023 


0.0227 ± 0.0004 






a c h 2 


0.11 


0.110 ±0.004 


0.11 


0.111 ±0.004 








0.75 


0.76 ± 0.03 


0.73 


0.74 ± 0.02 








0.003 


0.007 ± 0.005 










CL B 
"m 


0.25 


0.23 ± 0.03 


0.3 


0.3 ± 0.2 


0.1 


0.3 ± 0.2 


1Z A 


0.75 


0.76 ± 0.03 


0.4 


0.3 ± 0.2 


0.4 


0.3 ± 0.2 


X 


0.5 


0.3 ± 0.3 


0.7 


0.4 ± 0.3 


0.2 


0.5 ± 0.3 


(cj/H) 


6x 10" 10 


(3 ± 3) x 10" 10 


8 x 10" 10 


(4 ± 3) x 10" 10 


4x 10" 10 


(4 ± 3) x 10" 10 


a 


55.6° 


179.0° ± 101.0° 


37.0° 


162.8° ± 104.3° 


244.0° 


150.2° ± 102.5° 


P 


59.3° 


90.6° ± 38.3° 


26.4° 


82.1° ±39.6° 


106.9° 


75.3° ± 39.0° 


7 


102.6° 


180.5° ± 100.3° 


321.0° 


193.2° ± 103.8° 


170.3° 


203.8° ± 100.3° 
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(a) Flat-decoupled-Bianchi model 



-60.0 60.0 juK 

(b) Fixed-decoupled-Bianchi model 



Figure 10. Best-fit non-physical Bianchi VII/j templates found in full-sky ILC data. 
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